Formation of Liesegang patterns: Simulations using a kinetic Ising model 
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A kinetic Ising model description of Liesegang phenomena is studied using Monte Carlo simula- 
tions. The model takes into account thermal fluctuations, contains noise in the chemical reactions, 
and its control parameters are experimentally accessible. We find that noisy, irregular precipitation 
takes place in dimension d = 2 while, depending on the values of the control parameters, either 
irregular patterns or precipitation bands satisfying the regular spacing law emerge in d = 3. 
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I. INTRODUCTION 

Quasiperiodic precipitation patterns emerging in the 
wake of chemical reaction fronts are called Liesegang pat- 
terns |IJQ] . They have been studied for more than a cen- 
tury and a number of theoretical approaches have been 
developed to explain the experimental observations [||. 
Nevertheless, it was only recently that a model with input 
parameters fewer than the number of static and dynamic 
parameters characterizing the patterns has appeared [^ . 

This last theory is based on the assumption that the 
main ingredients of a macroscopic description should be 
a moving reaction front and the phase separation of the 
reaction product behind the front. Taking the proper- 
ties of the reaction front from the theory of the fronts 
in the A-\- B ^ C process and describing the phase 
separation process by the Cahn-Hilliard equation [^|, one 
arrives at a model with a minimal number of parameters. 

The above theory successfully explains that the po- 
sitions, Xn of the precipitation bands form a geometric 
series, a;„ ~ (1 -f p)" (spacing law 0), and gives the 
spacing coefficient p in terms of the initial concentra- 
tions of the reactants A and B in agreement with the 
Matalon-Packter law PJ|]. Furthermore, the parameters 
in the model can be determined from experiments and 
the time-scale of the emergence of a band can be calcu- 
lated [Q. Finally, the width law relating the position 
and the width of the bands can also be derived il^,|l^ in 
agreement with observations Q. 

The success and versatility notwithstanding, this the- 
ory needs further developments since, in its present form 
[Q , it is a mean field theory without the fluctuations be- 
ing accounted for. There are two ways to include the 
fluctuations. One is to add conserved thermal noise to 
the Cahn-Hilliard equation as it is done in Model B of 
critical dynamics |13| . In this case, there would be an 
additional problem of handling the noise in the reaction 
zone. It has been shown |1J], however, that this noise is 
irrelevant for d > 2 so, in principle, it could be neglected. 

In this work, another approach is taken for including 



the fluctuations. Namely, we shall study the kinetic Ising 
model version of the process and thus include the noise 
through the probabilistic description of the transitions 
between discrete states of the system. 

There are several reasons for our choice. First, the 
problem of difference between the handling of fluctua- 
tions in the diffusive and reaction processes does not 
arise. Second, the kinetic Ising version of the model has 
a meaning of a mesoscopic (or perhaps a microscopic) de- 
scription. Since the mechanism of band formation may 
be at work at a length-scale that makes possible the con- 
struction of submicron Liesegang structures, this kinetic 
Ising model approach may have a direct bearing on future 
experiments |T3]. Third, our choice was also influenced 
by having more expertise in simulations of kinetic Ising 
models. 

We shall start (Sec. II.) by a detailed discussion of the 
kinetic Ising model designed to describe the band forma- 
tion. The simulation results for this model are presented 
in Sec. III. First, the d — 2 case is treated where we do not 
find regular band formation. Then the d — Z simulations 
are discussed which show the emergence of Liesegang 
patterns satisfying the usual spacing law. Final com- 
ments about the difficulties of predicting temperature- 
dependent effects in Liesegang phenomena can be found 
in Sec.IV. 



II. KINETIC ISING MODEL DESCRIPTION 

A. General aspects of the theory 

The aim of the theories of Liesegang phenomena is to 
explain how a high-concentration electrolyte A diffuses 
into a gel soaked by a low-concentration electrolyte B 
and how the spatial distribution of the precipitate D is 
formed in the wake of the diffusive reaction front. Ac- 
cordingly, all the theories follow the scheme 



A + B^C 



D . 



(1) 



where C is an intermediate reaction product which is gen- 
erally not very well known. This uncertainty is then the 
basis for the existence of a number of competing theories 
||l6|-p3| with the differences arising from the interpreta- 
tion of C and from the level of details in the description of 
the dynamics of C-s. A significant drawback of all these 
theories is that they contain a large number of parame- 
ters and many of them are uncontrollable experimentally. 
Thus it is not entirely surprising that thorough compar- 
isons between experiments and theories have not been 
carried out. 



B. Cahn-Hilliard equation with source 

In our view, the uncertainty about the intermediate 
product and its dynamics can be used to advantage in 
building a general theory of band formation. One can 
interpret the concentration of C-s, c, as a kind of or- 
der parameter that takes a value Cp in the ordered (pre- 
cipitate) phase and another value Cg in the disordered 
(low-density) phase. The C-s are obtained from the 
A + B ^ C process so they are produced in the reaction 
zone. Furthermore, the dynamics of C-s obeys global 
conservation and it should be a phase-separation type 
dynamics since, in the expected final state, one has re- 
gions of high- (precipitation) and low-density (interband) 
regions in equilibrium. This phase-separation dynamics 
can be described on a coarse-grained level by the Cahn- 
Hilliard equation ||6| with the generation of C appearing 
as an additional source term. The resulting equation for 
the space- and time-dependent order-parameter density, 
c(x, t) is given by [0 



dtc = -AA 



Sfjc) 
6c 



tAc 



+ S{x,t). 



(2) 



Here A is a kinetic coefficient, /(c) is the Landau- 
Ginzburg free energy of the system which should have 
two equal minima at c = Cp and c = Cg. The term crAc 
with cr > provides stability against short-wavelength 
fluctuations, and finally S{x^ t) is the production rate of 
C-s in the reaction-diffusion process A + B ^ C. The 
properties of the source are known P,e3J23] . It is local- 



ized, its center, a:/, moves diffusively {xf = yJ^Dfi), and 
it leaves behind a uniform concentration cq of C-s. 

The parameters A and a in eq.(y) can be used to 
set the time-scale and length-scale, respectively, and the 
source, S", is completely specified by the initial densi- 
ties (ao,6o) and diffusion constants iJJ = Da ~ Db) 
of A and B j|j2j,|^. Only the function /(c) remains 
to be parametrized. One expects that the details of 
this function will not affect the overall properties of the 
pattern-formation process, the existence of two minima 
at Cp and Cg being the only important feature. Thus, as- 
suming e.g. a coexistence curve that is symmetric about 
c — {cp + Cg)/2, one can parametrize this function as 



/(c) = — e(c — c)^-)-7(c— c)^. The scale of the concentra- 
tion can be set by the parameter 7 and e remains a free 
parameter in the theory. Consequently, one has a theory 
which has only a single adjustable parameter apart from 
the parameters setting the scale of the length, time, and 
the concentration field. 

As discussed in the Introduction, all the observed gen- 
eral features of Liesegang phenomena can be derived from 
the above theory M , including the time-scale of the emer- 
gence of a single band and the length-scale for the width 
of the bands |13,p| . Actually, it is somewhat surprising 
that the mean- field level description in terms of eq.(||) 
performs so well. The reason for this may lie in the ex- 
perimental observation that the patterns are frozen (they 
do not evolve over time-scales extending up to 30 years 
) meaning that the dynamics takes place at a very low 
effective temperature, i.e. the noise is negligible. 

The noise may indeed be negligible in the late stages 
of the formation of precipitation bands but the initial 
stages should be related to some instabilities and there 
the fluctuations should play a more prominent role. In 
particular, the mean-field description relies on a spinodal- 
decomposition instability (the moving front generates 
particles and pushes the concentration past the spinodal) 
and the possibility for the precipitation to take place 
through a nucleation-and-growth mechanism (where the 
fluctuations are important) is completely lost. 



C. Kinetic Ising model with Kawasaki + Glauber 
dynamics. 



In order to include fluctuations in the Liesegang pro- 
cess, a kinetic Ising model will now be considered that is, 
we believe, a finite-temperature extension of the theory 
embodied in eq.(0). The model introduced below can be 
viewed in two ways. Either it is a discretization scheme to 
equation (0) and then the description is on a mesoscopic 
level with c{x, t) being the discretized concentration of 
the order parameter, C. Or, it can also be viewed as a 
simulation of the stochastic motion of C-s which are now 
particles at microscopic scales (in this case the coarse 
graining has been carried out in time). In the following 
we shall use the latter "particle" language. 

Let us begin the introduction of discretized description 
by identifying the particles C with the up-spins of an 
Ising model on a hypercubic lattice. Then the down-spin 
sites are empty places and the formation of precipita- 
tion bands is modeled by a combination of spin-flip and 
spin-exchange dynamics P6| . Namely, the initial state 
is prepared with all spins down (no C particles present) 
and the localized front flips the down spins (Glauber dy- 
namics p7|) thus producing the C-s. The flip rate Wf at 
site r is given by 



wp^ S{x,t)-{1 - ap) , 



(3) 



where a^ = ±1 is the Ising spin at site r = {x,f±) with 
X being the coordinate in the direction of the motion of 
the reaction front while r± representing the coordinates 
in the transverse direction (length is measured in units 
of the lattice spacing a). The factor (1 — ap) ensures that 
the front flips only down spins (the particles are pro- 
duced in the front and the back reaction is negligible). 
Finally, S{x,t) is a function describing the motion of the 
reaction front and the change of the reaction rate with 
time. The front is assumed to be homogeneous in the 
transverse direction and its actual shape, S{x,t), can be 
taken from the solution and simulations oi A + B ^ C 
process PJ28|. Since the width of the reaction zone is 
small and it changes with time very slowly, for all prac- 
tical purposes, the function S{x,t) can be approximated 
by a constant within a small interval [xf — A,Xf + A] 
around the center of the reaction zone Xf — ^JID ft 

A 

S{x,t)^^{x-Xf + A)0{xf + A~x). (4) 

where 9{x) is the step function and the amplitude 



A V^ 



(5) 



is chosen such that the front leaves behind a constant 
(co) concentration of particles ^jg. 

Once the particles are created, they diffuse and inter- 
act. This part of the dynamics can be described by a 
spin exchange process (Kawasaki dynamics pS]). The 
rates of the exchanges are assumed to satisfy detailed 
balance at temperature T with ferromagnetic coupling 
(J > 0) between the spins in order to describe the ex- 
pected attraction among the C"s. Assuming the usual 
nearest-neighbor Ising Hamiltonian 



H^-J 






(6) 



the rate of exchange between neighboring sites f and r ' 
can be chosen to be IM 



1 



1 + gSE/iksT) 



(7) 



where Tg sets the timescale, T and fc^ arc the tempera- 
ture and the Boltzmann constant, respectively, and 6E 
is the change in the energy 6E — 5H{(jp ^^ ap) due to 
the exchange of spins at r and r '. 

Without the spin-flip dynamics, the system would re- 
lax to the equilibrium of the Ising model at temperature 
T and at fixed magnetization. Thus choosing the tem- 
perature low enough (below Tc of the Ising model) and 
making the spin-flip front produce the right magnetiza- 
tion density, the system will be in the unstable part of 
the phase diagram of the Ising model and phase separa- 
tion will take place. If this model represents the pattern 
forming process correctly then one expects the emergence 
of bands of up and down spins in the wake of the moving 
spin-flip front. 



III. SIMULATION RESULTS 



A. Parameters 

Monte-Carlo simulations of the kinetic Ising model de- 
scribed above have been performed in dimensions d=2 
and 3. The dimension-specific properties of the lattices 
used in the simulations will be given in the appropri- 
ate subsections below. Here we enumerate and discuss 
only those adjustable parameters which are used inde- 
pendently of dimension: 

• The temperature T is measured in units of J/ks, 
where J is the nearest-neighbor coupling of the 
Ising model (||). It is clear from the considerations 
of the previous section that T < Tc should be used. 

• The particle concentration, cq, deposited by the 
front must be chosen so that, at the given T, it 
places the system in the metastable or unstable re- 
gion of the phase diagram of the Ising model. 

• Length is measured in units of the lattice spacing a. 
It should be noted that the value of a depends on 
the interpretation of C and a can be a microscopic 
length-scale as well as a mesoscopic one. 

• Time is measured in units of the "microscopic" 
timescale Tg [see eq.(0)]. Again, this timescale may 
be coming from microscopic or mesoscopic pro- 
cesses depending on the interpretation of C. 

• The diffusion coefficient of the front (Df) is, in 
principle, a well defined quantity [H|. The uncer- 
tainty of the connection between the microscopic 
and macroscopic length- and time-scales, however, 
makes it difficult to set a value for Df in terms of 
a and t^. We shall thus treat Df as a parameter 
that can be freely varied. 

• The width of the reaction front (2A) does not ap- 
pear to influence the emerging patterns (this is the 
experience both from the mean-field theory jj] and 
from small scale simulations). Thus, in most of our 
simulations, we set A — 1/2 i.e. the front coincides 
with one of the lattice planes perpendicular to the 
motion of the front. 

The three important parameters are T, co and Df, and 
one should search for pattern formation in this three di- 
mensional parameter space. This is not necessarily an 
easy task since it is known, both experimentally and 
theoretically ||2^, that patterns are formed only in a re- 
stricted domain of the available free parameters. Since 
our search is finite, our statements (especially about the 
absence of patterns) are always pertinent only to the pa- 
rameter domain investigated. 



B. Two-dimensional simulations 



Simulations have been performed on stripes of length 
Lx and width Ly = 50 or 100 with periodic boundary 
conditions in the transverse (y) direction. L^ was chosen 
to be long enough so that no end effects would be ob- 
served with Lx = 3000 being a typical value. The initial 
position of the front was always at a; = 1 and we used 
closed boundary conditions (no particle crossing) at both 
X — and x = L^- 

The possible formation of patterns was investigated 
for the following ranges of the adjustable parameters: 
0.02 < T < 1.5, 0.01 < Df < 1 and 0.05 < co < 0.5. 
About 200 sets of values have been studied in the above 
domain and no regular patterns were observed. A typical 
result is displayed in Figjl] for T = 0.7, £>/ = 0.025 and 
Co = 0.3. 



Thus the Cahn-Hilliard description [eq.(H)] should apply 
and we should observe regular band formation. 

The nonexistence of regular banding should not be con- 
sidered as a contradiction with the experimental obser- 
vations oi d = 2 Liesegang patterns. The experimental 
systems always have macroscopic width in the third di- 
mension and this appears to stabilize the patterns. In- 
deed, Fig.^ shows a simulation with the same parameters 
as those in Fig.y, except an extra layer in the third dimen- 
sion is added. As one can see, the bands in the two-layer 
system are much better defined and they display some 
regularity. Note in particular, that the concentration 
within the bands in Fig.0 has reached the equilibrium 
value (c ~ Cp ~ 1) while the maximum concentration 
regions in Fig.n^ are roughly halfway in between the equi- 
librium values, c„ « and Cp ~ 1. 
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FIG. 1. Concentration profile of C-s averaged over the 
transverse direction for a two-dimensional pattern obtained 
in a stripe of size 50 x 3000 with T = 0.7, Df = 0.025 and 
Co = 0.3. The distance along the slab (a;) is measured in units 
of the lattice spacing. The dotted line shows the position of 
the front at the time (i = 1.8 * 10 ) the concentration was 
measured. The time is in units of the inverse of the rate of 
hopping for free C particles (r^). 



FIG. 2. Concentration pattern in a system of size 
50 X 2 X 3000. The parameters T, Df and Cq and the no- 
tations are the same as on Fig.nl. The dotted line shows the 
position of the front. 



C. Three-dimensional simulations 



We have also made a few simulations outside of the 
above domain in order to check the possible strong ef- 
fect of the changes in a single parameter. These nonsys- 
tematic searches did not lead to pattern-forming regimes 
either. 

We have extended the simulations to cases where the 
interactions are not restricted to nearest neighbors but 
extend up to seven lattice spacings. No patterns were 
found, although the randomness of the pattern slightly 
decreased when the range of the interaction was in- 
creased. This is expected as in the limit of long-range in- 
teractions one should reach a continuum, noiseless limit. 



The simulations were performed on slabs of length Lx 
and of cross-section of size Ly x L^- Periodic boundary 
conditions were used in the transverse (y, z) directions 
while closed boundary conditions were employed at the 
two ends of the slabs. As in the two-dimensional case, L^ 
was chosen so {L^ ~ 3000) as to avoid end effects from 
the transverse wall a,t x = L.^- Most of the simulations 
were done for Ly — L^ = L — 10, 20 and 40 in order to 
observe finite-size effects. The initial condition was again 
an empty (all spins down) state with the front situated 
at a; = 1. 
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FIG. 3. Time evolution of the concentration of C particles 
averaged over the transverse {y, z) directions in slabs of size 
20 X 20 X 3000. The parameters (T = 1.0, D/ = 1.0, co = 0.4) 
were chosen to be in the "coarsening pattern" regime. The 
time is given in units of Te that is the inverse of the rate of 
hopping for free C particles. The distance along the slab (x) 
is measured in units of the lattice spacing. The dotted lines 
show the positions of the reaction front at time t. 



Precipitation patterns were observed for T < 1 if 
cg > 0.25 was chosen well in the metastable or unsta- 
ble region of the phase diagram of the d=3 Ising model. 
For large front diffusion [Df > 1), these patterns were 
not regularly spaced and were not stable. Namely, coars- 
ening was observed within reasonable observation time 
(time of formation of about 7-10 bands). A typical ex- 
ample is shown on Fig|3[ 

The bands are more stable at lower temperatures (T < 
0.8) and their spacing becomes more regular, as the front 
diffusion coefRcient is decreased below Df < 0.1. For 
Df < 0.05, one observes the emergence of Liesegang type 
band patterns, typical examples being those shown on 
Fig.^-c. In this case, the system is of size {Lx Lx 3000) 
and the parameter values used are T = 0.7, Df = 0.025 
and Co = 0.3. Results of runs for three cross sections 
{L = 10, 20 and 40) are displayed. Comparing these pic- 
tures, one cannot see any obvious finite-size trends. 

In order to investigate the spacing law one would need 
a large number of bands. Unfortunately, in the regime 
where the best Liesegang type patterns are obtained one 
is restricted in the extent of explorations by the com- 
puting resources. Due to the low value of Df, the front 
is moving very slowly and consequently the computation 
time for obtaining e.g. 10 bands becomes very large. The 
CPU time necessary to produce the pattern on Fig.^ 
was about 1500 hours on a Sun Ultra- 10 workstation. 

Fig.H shows roughly the limits of the possibilities of 
our simulations at present. To establish the spacing law 
firmly, we would certainly need more bands. Using the 
last 4-6 bands obtained from pictures similar to those on 
Fig.H, one can see that the positions of the bands (xn) do 
approximate a geometric series Xn ^ (1 +p)" well, and 
one can extract an approximate the spacing parameter, 
p. In general, we find that the spacing coefficient does 
not show discernible finite-size trends. For example, the 
spacing parameter is p ~ 0.17 for all values of L in Figjj. 

An interesting feature of the patterns emerging in our 
simulations is the presence of material (C) in between 
the bands (see Figures pla-c and Ha-c). Of course, it is 
not entirely clear whether part of this material should 
be considered as a low density precipitate (seen in many 
experiments ||]) or should it be just regarded as a "gas" 
phase of particles C. Visual inspection of the interband 
region reveals the presence of both small clusters and sin- 
gle particles. The small clusters live long (especially at 
lower temperatures), their lifetime is comparable to the 
time of formation of several bands. Thus the interpre- 
tation of part of the interband material as low density 
precipitate may have some validity. In order to make 
progress in this problem, one would need larger scale 
simulations as well as more understanding of the con- 
nection between the microscopic and macroscopic time- 
and length-scales. 



IV. FINAL REMARKS 
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Apart from the simplicity, an important feature of our 
model is that fluctuations are included. A frequent conse- 
quence of the presence of fluctuations is the disapperance 
of order in low dimensions and, indeed, our simulations 
also indicate that Liesegang type patterns are absent in 
d = 2 dimension while they exist in d = 3 dimensional 
samples. 

Another advantage of including the noise is that the 
model has pattern-forming regimes which have charac- 
teristics of cither the prenuclcation theories [po|,pi[ or 
the postnucleation competitive growth theories ]22|. In- 
deed, one can see (Fig.Qb) that the bands are forming 
at the position of the reaction front for Df small, while 
the bands are formed as a result of coarsening and com- 
petitive growth well behind the front in case of large Df 
(Fig.g) . Note that the absence of Liesegang type patterns 
in the second regime is in agreement with the inability 
of producing such patterns in postnuclation competitive 
growth theories 22 1. 

We feel that the most important feature of the model 
is that it makes clear that the important and experimen- 
tally controllable parameters are T, cq and Df. Indeed, 
both Co and Df are known functions of the initial den- 
sities (ao,6o) and diffusion constants {D = Da ~ Df,) of 
reagents A and B in the A + B ^ C reaction P,e3,E5[ . 
Although the diffusion constants Da and Df, are usually 
not controllable, qq and 6o can be set to given values and, 
consequently, cq and Df can be varied independently. 

Among the three parameters, a change in T leads to 
unpredictable changes in the various diffusion coefficients 
and background processes in a real Liesegang experiment 
(for example, the diffusion coefficients of the background 
ions may also be important |Q in determining the spac- 
ing coefficient). Thus it is advisable to keep the tempera- 
ture constant, and the most promising and, indeed, most 
often used way of looking for trends in experiments PJ9] 
is the change of the concentrations of the inner and outer 
electrolytes (ao and bo). 

An obvious consequence of our simulations is that one 
should be able to observe a crossover from prenuclcation 
regime to a regime where postnucleation processes are 
dominant. Namely this could be achieved by tuning ao 
and bo in such a way that, at fixed co, the front diffusion 
coefficient Df would be varied over as much a range as 
possible. 



FIG. 4. Concentration averaged over the transverse {y, z) 
directions in slabs of size L x L x 3000. The values of the 
parameters (T — 0.7, Df = 0.025 and co = 0.3) are the same 
for three transverse sizes L = 40 (a), L = 20 (b) and L = 10 
(c). The distance along the slab (x) is measured in units of 
the lattice spacing. The dotted lines show the positions of the 
reaction front (the front is off the scale on Fig.Mc). 



ACKNOWLEDGMENTS 

This work has been supported by the Swiss National 
Science Foundation and by the Hungarian Academy of 
Sciences (Grant OTKA T 029792). 



[1] 

[2] 
[3] 



[4] 

[5] 
[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[12] 



[13] 

[14] 
[15] 



[16] 
[17] 



R.E. Liesegang, Naturwiss. Wochenschr. 11, 353 (1896); 
H. K. Henisch, '^Periodic precipitation" , Pergamon Press, 
(1991). [18] 

For a recent comparison between theories and experi- 
ments see T. Antal, M. Droz, J. Magnin, Z. Racz and [19] 
M. Zrinyi, J. Chem. Phys. 109, 9479-9486, (1998). [20] 

T. Antal, M. Droz, J. Magnin, and Z. Racz, Phys. Rev. [21] 
Lett. 83, 2880 (1999). 

L. Galfi and Z. Racz, Phys.Rev.A38, 3151 (1988). [22] 

J. W. Cahn and J. E. Hilhard, J. Chem. Phys. 28, 258 
(1958); J. W. Cahn, Acta Metall, 9, 795 (1961). 
K. Jablczynski, BuU. Soc. Chim. France 33, 1592 (1923). 
R. Matalon and A. Packter, J.CoUoid Sci. 10, 46 (1955); 
A.Packter, KoUoid Zeitschrift 142, 109 (1955). 
Z. Racz, Physica A274, 50 (1999). 
J. Magnin, PhD Thesis, University of Geneva, 2000. 
S. C. MiiUer, S. Kai, and J. Ross, J. Phys. Chem., 86, 
4078 (1982); B. Chopard, M. Droz, J. Magnin and M. [23] 
Zrinyi, J. Chem. Phys. 110, 9618 (1999). 
P.C. Hohenberg and B.I. Halperin, Rev. Mod. Phys. 49, 
435 (1977). 

S. Cornell and M. Droz, Phys. Rev. Lett. 70, 3824, (1993) 
Up to few years ago, the smallest distances between [24] 
bands achieved in experiments were of the order of A ~ 
10/i. Recent reports show a decrease of A to A « 1/^ [X.G. [25] 
Zheng, M. Taira, M. Suzuki, and C.N. Xu, ApL Phys. 
Lett. 72, 1155 (1998); P. Hantz, Pattern Formation in 
Reaction- Diffusion Systems (Thesis, Eotvos University, [26] 
Budapest, 1998).]. 

W. Ostwald, '^Lehrbuch der allgemeinen Chemie" , Engel- [27] 
man edt., Leipzig (1897). [28] 

C. Wagner, J. Colloid Sci., bf 5, 85 (1950); S. Prager, 
J. Chem. Phys. 25, 279 (1956); Ya .B. Zeldovitch, G.L [29] 
Barrenblatt and R.L. Salganik, Sov. Phys. DokL 6, 869, [30] 
(1962). In these theories the reactants turn directly into 



precipitate {A + B ^ D) thus the dynamics of the inter- 
mediate product C does not arise. 

N.R. Dhar and A.C. Chatterji, KoUoid-Zeitschrift 37, 2-9 
(1925); J. Phys. Chem. 28, 41-50 (1924). 
S. Shinohara, J. Phys. Soc. Japan 29 1073 (1970). 
G.T. Dee, Phys. Rev. Lett. 57, 275 (1986). 
B. Chopard, P. Luthi and M. Droz, Phys. Rev. Lett. 72, 
1384 (1994); Jour. Stat. Phys. 76, 661-676 (1994). 
M. Flicker and J. Ross, J. Chem. Phys., 60, 3458-3465 
(1974); S. Kai, S. C. MiiUer, and J. Ross, J. Chem. Phys., 
76, 1392 (1982); R. Feeney, S.L. Schmidt, P.Strickholm, 
J. Chadam, and P. Ortoleva, J. Chem. Phys., 78, 1293 
(1983); G. Venzl, J. Chem. Phys., 85 1996, 2006 (1986). 
These works attribute the Liesegang band formation to 
the instability present in the post-nucleation coarsening 
process. They cannot explain, however, the regular as- 
pects of the band spacing. 

M. Chacron and I. L'Heureux, Phys. Lett. A 263, 70 
(1999); I. L'Heureux, Phys.Rev.E62, 3234 (2000). This 
is an attempt to make a synthesis of pre- and post- 
nucleation theories. No quantitative analysis of the re- 
sulting patterns were carried out. 

The case of unequal diffusion coefficients is treated in Z. 
Koza, J. Stat. Phys, 85, 179 (1996). 

More complicated reaction schemes {A + 2B -^ C and 
2A + B ^C) are studied in J. Magnin, Eur. Phys. J. B 
to appear. 

M. Droz, Z. Racz, and J. Schmidt, Phys.Rev.A39, 2141 
(1989). 

R.J. Glauber, J. Math. Phys. 4, 294 (1963). 
H. Larralde, M. Araujo, S. Havlin and H. E. Stanley, 
Phys. Rev. A46, 855 (1992). 
K. Kawasaki, Phys.Rev.145, 224 (1966). 
T. Unger and Z. Racz, Phys.Rev.E61, 3883 (2000). 



